The Conal representation of Quantum States 
and Non Trace-Preserving Quantum Operations 
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We represent generalized density matrices of 2 a d-complex dimensional quantum system as a sub- 
cone of a real pointed cone of revolution in R d , or indeed a Minkowskian cone in E 1,d - 1 . Gener- 
alized pure states correspond to certain future-directed light-like vectors of E 1,d _1 . This extension 
of the Generalized Bloch Sphere enables us to cater for non-trace-preserving quantum operations, 
and in particluar to view the per-outcome effects of generalized measurements. We show that these 
consist of the product of an orthogonal transform about the axis of the cone of revolution and a 
positive real linear transform. We give detailed formulae for the one qubit case and express the 
post-measurement states in terms of the initial state vectors and measurement vectors. We apply 
these results in order to find the information gain versus disturbance tradeoff in the case of two 
equiprobable pure states. Thus we recover Fuchs and Peres' formula in an elegant manner. 
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I. INTRODUCTION 



The space of pure states of finite d-dimensional Quan- 
tum Mechanics CP d , set of rays in the complex Hilbert 
space C d , is, as most complex spaces, not easy to vi- 
sualize. Physical motions, let alone unitary time evolu- 
tions, have no clear geometric interpretation. However, 
the set of hermitian operators on C d , Herm^C), is a 
(d 2 — l)-dimensional real vector space, and as such is cer- 
tainly easier to represent geometrically. States, or more 
generally density matrices, form of course a subset of 
Herm,j(C). The Generalized Bloch Sphere Representa- 
tion (U|-|a]) is a famous application of this fact which 
has proved to be popular and elucidating: a given den- 
sity matrix can be represented as a real vector inside a 
(hyper) sphere. 

It turns out that this representation defined for density 
matrices, or unit trace positive elements of Herm^C), is 
only good at handling unitary or trace preserving quan- 
tum operations on density matrices: the former induce 
rotations of the Bloch vector, the latter affine transfor- 
mations 0]. Individual outcomes of generalized mea- 
surements, for example, are not directly representablc. 
Considering the insight the Bloch Sphere representation 
gave to unitary and trace-preserving operations, it seems 
interesting, for the mere sake of geometry at first, but 
mainly to give a useful picture to tackle Quantum In- 
formation problems, to extend it to cater for non-trace- 
preserving quantum operations. This is further moti- 
vated by the fact that the space of (semi-definite or def- 
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inite) positive hermitian operators, thereafter denoted 
Herm+(C), is a closed convex cone, and that all ad- 
missible quantum operations should be a subset of the 
transformations of this cone. 

In spite of being so central in Quantum Information 
Theory, the tradeoff between how much Shannon Infor- 
mation one may gain about a quantum system versus 
how much Disturbance the observation must necessar- 
ily cause to the system, remains extremely difficult to 
quantify. Quantum cryptographists tend to circumvene 
the problem: most of their proofs are a witty blend of 
the particular symmetries of the protocol in question, 
together with a convoluted machinery. A few attempts 
have been made to solve the tradeoff Q p| , but only one 
deals with discrete ensembles - namely the case of two 
equiprobable pure states, and this is already something. 
Unfortunately the approach involves lengthy algebra and 
a number of assumptions. One should be able to find a 
method which gives a glimpse of intuition about the ge- 
ometry of optimal measurements, and for this purpose, 
we think that our approach is useful. 

In section[n]we consider general quantum systems of d 
complex dimensions. We give a representation of the set 
of positive hermitian matrices HermJ (C) as a subcone of 
a real Minkowskian cone in R d2 , and analyse geometrical 
properties of generalized measurements in this setting. 
We find that our approach is particularly useful to rep- 
resent per-outcome post-measurement states, and that 
pure states correspond to certain light-like vectors of the 
cone. Unitary operators on the complex system become 
real orthogonal transforms, while positive operators be- 
come real positive transforms. Section II I II should be of 
special interest for quantum information theorists: we 
treat the d = 2 one qubit case in full detail. We find 
further geometrical relations between measurement vec- 
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tors, state vectors and post- measurement state vectors 
and give explicit formulae. In section II VI we apply our 
results to a typical Information gain versus Disturbance 
tradeoff scenario in which Alice gives Eve two equiprob- 
able pure states. Thus we recover Fuchs and Peres' for- 
mula in an elegant and geometrical manner. 



II. CONAL REPRESENTATION OF 
d-DIMENSIONAL QUANTUM SYSTEMS 

The state of such a system is described by a d x d 
density matrix. We shall express hermitian matrices as 
real linear combinations of Hilbert-Schmidt-orthogonal 
hermitian matrices, and then restrict this representation 
to elements of Hermj"(C), or generalized density matri- 
ces. Herm^(C) turns out to be "isomorphic" to a con- 
vex subcone of a cone of revolution in R rf , or indeed a 
Minkowskian future cone in E 1 ^ . We then analyse 
the effects of quantum operations on density matrices in 
this representation. 



A. Hermitian matrices 

Let {n}, i € {1, . . . , d 2 — 1}, be a Hilbert-Schmidt 
orthogonal basis (as in Q)) of d x d traceless hermitian 
matrices, and let To be the identity matrix I. Throughout 
this article latin indices will run from 1 to d 2 — 1, greek 
indices from to d 2 — 1, and repeated indices are summed 
unless specified. We take the r M 's to satisfy by definition: 

V>, v Tr(-r M T„) = dS^ (1) 

with 5 the Kronecker delta. {t m } m is a basis of Herm t ;(C), 
and any hermitian matrix A 6 Herm^C) decomposes on 
this basis as 

A= ~(Tr{A)I + Tr(ATi)Ti) 

= ~Tr(Ar M )T M (2) 

Letting A = Q4J S R d ' with A^ = Tt(At^) be the 
component vector of A in this particular basis, we have 

Vi,5e Herm d (C), AB = ^A^B v t^t v 

hence TrAB = -A.B = -,A n B n (3) 
d d~^~^ 

2 ^ 

We shall call A the vector in R d , A = (AA the restricted 
vector in R , and </> the coordinate map: 

4> : Hcrm d (C) -> R d ' 2 
A^A 



Equation says that <f> is an isometric isomorphism of 
(Hcrm d (C),Tr( )) onto (R d2 , (l/d)( . )). Therefore any 
linear operator L on Herm^C) defines via </> and <f>~ x an 
operator on R d , M(L) = cp o L o ^> _1 . This definition 
yields the following "morphism" property : 

Lemma 1 If Li,L2 are linear operators on Herrrid(C), 
then M(Li) = </>oLjO0 _1 for i= 1,2 are endomorphisms 

,2 

of M. and satisfy 

M(Lt o L 2 ) = M(ii)M(i 2 ) (4) 

In particular, any complex d x d matrix A defines via 
Ad a ■ p i— > ApA^ a linear operator on Hermd(C) which 
corresponds to a real endomorphism M(AdA) '■ p 
M(Ad A )p ; and Ad AB = Ad A oAd B implies M(AcUb) = 
M(AdA)M{AdB)- As a direct consequence of this and 
the previous definitions , calling GL n (K) the group of 
invertible n x n matrices on the field K, we get : 

Lemma 2 For any subgroup G of GL^C), the following 
mapping 

^:G^tp(G)cGL d 2{R) 

A i ► M(Ad A ) = o Ad A o (jy (5) 

is a group homomorphism. ipiG) is a subgroup of 
GL d 2 (R) . 

Note that since ^(1) = tjj(— I) = I, ijj is not necesseraly 
injective. Moreover tp is certainly not linear. An inter- 
esting subgroup is the Special Unitary group SU(d) = 
{U S GL d (C) I UW = I, detU = 1}. We call 
SO(n) = {O G GL n (R) I 00 T = I, det = 1} the 
special orthogonal group in n-dimensions. 

Lemma 3 Special Unitary transformations on 
Herm d (C), Ad v : p h-> UpW with U e SU(d), in- 

,2 

duce rotations of R about the l-axis. In fact the linear 
transforms ip(U) : p i— > ip(U)p are special orthogonal and 
ip(SU(d)) is a subgroup of SO(d 2 — 1). It is a proper 
subgroup when d>3. Moreover, ip(U(d)) — tjj(SU(d)). 

Proof: Let p = (l/d)(Tr(p)I + p^i) a hermitian matrix. 
Using © and the fact that Adu is trace-preserving for 
U unitary: 

upu f . u P u f = (Tvpf + {m)p)M u )p)i 

= dTi{UpU^U P U^) 
= dTip 2 = p.p 
= (Trp) 2 +P l P l 

In addition to preserving the first component p Q = Trp, 
ip(U) preserves the R d2_1 scalar product p. -p.- For all 
U G SU(d), there exist t g R and B e su(d) such 
that U = U(t) = exp(tB). Since det -0(^(0)) = 1 and 
t i ► det ip(U(t)) is continuous and has values in {±1}, 
detip(U) = 1. Thus ip(U) is a special rotation about the 
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I -axis of R d . By LemmaH if>(SU(d)) is a subgroup of 
the special orthogonal group SO(d 2 — 1) C SO(d 2 ). As 
for all 6 G R, Ady = Ad e<9[ /, ^(U(d)) = ip(SU(d)). 
Since the {\/— It,} span the Lie algebra su(d), ip(SU(d)) 
is the Adjoint group of SU(d). For d — 2, we get the 
whole of 50(3), but this is not the case for d > 2, as is 
easily seen looking at the dimensions: 

dim SU(d) = d 2 - 1 
dim SO(d 2 - 1) = i(d 2 - l)(d 2 - 2) 

and dim SU(d) < dim SO(d 2 - 1) for d > 2. □ 
All the results of this section remain true of course when 
we just consider HermJ(C). From now on, for any A 
complex d x d matrix we shall denote ip(A) — M(Ad^) 
the real endomorphism of K d . 



B. Generalized density matrices 

In 2] , Zanardi showed using a restriction of a mapping 
analogous to 4> '■ A i— ► A that d x d density matrices lie 
m a convex subset of a ball ScR" . We shall extend 
this to a convex cone by considering generalized density 
matrices, by which we mean elements of Hermj~(C). This 
seems more natural in the sense that we like to think of 
the space of states of most physical theories and indeed 
Quantum Mechanics as a space invariant under positive 
linear combinations and not just convex combinations. 
In addition, this bigger space allows a per-outcome rep- 
resentation of generalized measurements. 
We define generalized pure states to be generalized den- 
sity matrices which yield pure states after rescaling them 
to unit trace. Note that these are not the "states of par- 
tial purity" of the complex d-dimensional system, which 
are singular density matrices. In other words, general- 
ized pure states are not the elements of the boundary 
of Herm^(C) in the sense of characteristic functions of 
cones (see for example). 

Proposition 1 The cone of positive hermitian matrices 
Herrri^ (C) is isomorphic to a convex subcone C of the 

,2 

following cone of revolution in K : 
d 2 -i 

T = {(A M ) G M d V Yl A < < ( d - !) A o, A o > 0} (6) 

i=l 

The set of generalized pure states verifies C = C PI dY , 
where dY stands for the boundary ofT. 

Proof : We begin as in 0. Let f denote the space of 
(not generalized) pure states in Hemic; (C). In addition 
to being positive, p € ? satisfies Tr(p 2 ) = Tr(p) = 1, so 



we have 

= j( 1+ £!i£i) = hence 
EiEi = d ~ 1 (7) 

The restricted vector (p.) is on a (d 2 — 2)-sphere of ra- 
dius \/d — 1, dS d ~ 2 , where S is the corresponding ball. 
In M. d , p pure sits in the intersection of the cylinder (JJJ 
and the p Q = Tr(p) = 1 hyperplane, in other words on 

dS d2 - 2 "centered" at (1,0,..., 0). 

Any density matrix can be expressed as a positive (con- 
vex) linear combination of pure states, and any posi- 
tive (convex) linear combination of pure states defines 
a density matrix. Calling D the set of (not general- 
ized) density matrices, D C Hull(J') and Hull(T) C D. 
Since D is closed, D = HuU(J'), a well-known result. 
As 4> : Herm c i(C) — * R d is linear and bi-continuous, 
<j>{D) = 0(Hull(V)) = 0(Hull(T)) = Hull(0(T)). This set 
is a closed convex subset of S "centered" at (1, 0, . . . , 0): 

<f){y) c dS d2 - 2 ^ HuU(^(T)) c 3 

Calling S + = Hull(0(CP)) the image set of density matri- 
ces as a subset of R d -1 , we get: 

pe^^Tr(p)=p = land (p.) € 5+ 

Now a non-zero p € Herm c ;(C) is positive if and only if 
(1/Trp)p is positive, that is if and only if ((1/Tr(p))p ) € 
S + . In R d2 , recalling that Tr(p) = p Q , this reads 

p G Herm+(C) <=>/)£ {(A , (A,)) G M d V(A,) G A 5+} 

(8) 

This clearly defines a cone C in ]R d2 . As S + C S 1 , C is 
a subcone of the cone of revolution T given by ©. 4> 
being an isomorphism, C is convex and isomorphic to 
HermJ(C). As pure states correspond to some points 

on the sphere dS d ~ 2 , generalized pure states lie in the 
boundary of V. Calling C the set of vectors of C corre- 
sponding to generalized pure states, we have C C CndT. 
Moreover C D C n dT follows from the fact that any 
rescaled positive matrix p such that Tr(p 2 ) = Trp = 1 
is a pure state. Remember that C is not the boundary 
of C, but a cone over 4>(7), the image set of pure states. □ 

As we shall see in detail in section Mil in d = 2 di- 
mensions, generalized pure states correspond to future- 
directed light-like vectors of Minkowski space of signa- 
ture (1,3). We have shown that this remains true to a 
certain extent in d-complex dimensions, Y being the fu- 
ture light-cone of Minkowski space E 1,d _1 with metric 
i]pv = Diag(d — 1, — 1, . . . , —1). Thus the appearance of 
a Minkowski product is to be expected. 
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As a consequence of Lemma |3 unitary transforms, 
since they leave HermJ(C) invariant, yield rotations 
which leave C (globally) invariant. This fact deserves 
to be analysed in detail to understand the geometry of 
C . As Unitary transforms act transitively on pure states, 
S + is the closed convex hull of a homogeneous subspace 
0(CP) of dS d ~ 2 . For the moment however, we shall con- 
sider the geometric representation of general quantum 
operations in C . 



C. Generalized measurements 

We call a generalized measurement Q a finite set 
{M m } m of complex d x d matrices which satisfy : 
E™ MlM m = I The set of {E m } m = {MlM rn } rn de- 
fines a Positive Operator Valued Measure (POVM), as 
E m € HermJ(C) and ^2 m E m = I. Given a quantum 
state or density matrix p 6 D, the generalized measure- 
ment {M m } m on p yields outcome m with probability 
p(m) — Tr(E m p), and if outcome m occurs, the post 
measurement state is p' m = (l/Tr(£' m ,o))(M m pM I | i ). We 
shall call p m = M m pM m £ HermJ(C) the unrescaled 
post-measurement state. 

Recall that any complex matrix can be polar-decomposed 
into a product of a unitary matrix and a positive matrix. 
For all m, there exists U m £ U(d) and A m £ HermJ(C) 
such that M m = U m A m . As E m = M m M m = A m A m , 
A m = v 'E m , the positive square root of E m . Using this 
polar decomposition, p' m is represented in the cone C by 

Pm = </>(Pm) = T^jk T(t>{U m ^E~p^jE\nU m ) 

Tr(£ m p) 

Tr (V E m py/ E m ) 

Thus when outcome m occurs, the post-measurement 
state p' m of {M m } m is the same as that of {y/E m } m up 
to a rotation ip(U m ), and similarly for the unrescaled 
states. As a consequence we shall consider the geometri- 
cal effects of generalized measurements {\JE m } m where 
E m and \JE m are in Herm+(C) and £ m E ™ = I bear- 
ing in mind that the most general measurements just 
involve rotations on the post-measurement state vectors. 
For example, in section IIVI Eve is free to perform uni- 
tary transforms on her post-measurement states, and can 
decide this according to the outcome m. The procedure 
we use to find the Disturbance is to first measure with 
{\fEm\m and then maximise on Unitary transforms act- 
ing upon post-measurement states. Using the conal rep- 
resentation, both sets of vectors { E m } m and {>/E m }m 
are in C, and ^2 m Em — (d, 0, . . . ,0). This enables us 
to represent elements of a measurement inside C, and 
visualize the action of a particular non-trace-preserving 
operation \JE m on a given density matrix p, in other 
words find p m in terms of E m or y/E m . 



D. Quantum operations represented in C 

One might wonder here why not just rescale all the 
post-measurement states and only consider the density 
matrices p' m . The reason for not doing so is that the un- 
rescaled states encode extra information: their "height" 
in the cone, the first component p mQ — Tr(E rn p), is sim- 
ply the probability of their outcomes. Under a given gen- 
eralized measurement, post-measurement vectors with 
identical first components are equiprobable. Thus the 
sections of C of constant Ao have a clear physical inter- 
pretation. We shall need the following simple properties: 

Lemma 4 For A £ Hermd(C) and B,C £ Herm^iC), 

Tr(BC) > 
Tr(BABA) > 

Proof: Let B = VB\Tb, then Tr(BC) = 
Tr^/BCVB) > since \fBCj~B £ Herm+(C). Then 
polar decompose A into A = U\A\, with U unitary and 
\A\ £ Herm+(C). As A £ Herm d (C), A = \A\ W '= A\ 
and 

Tt(BABA) = Ti(BU \A\ B \A\ ?7 t ) = Tr(U*BU \A\ B \A\) 

This is non-negative by the previous result since 
U^BU,\A\B\A\ GHermj(C). □ 

Unitary transforms induce rotations in C, and gener- 
alized measurements have the following geometric prop- 
erties: 

Proposition 2 The linear transforms ip(\/E m ) : p i— > 
p m associated to a generalized measurement {y/E m } m 
correspond to real symmetric matrices which are positive. 
They individually map C into itself. In addition, for any 
generalized pure state 9, ip(6) maps C into C. 
The probability of outcome m for a quantum system in 
state p is given by 

p{m) = 2 Em - p ( 9 ) 
Proof : By using @ successively, we have 

Pm = ^{\f E ~niP^f E ~mTp) 

= -Tr(v / ^T I/ v / £^7>)p i/ (10) 
^ M ^ p v 

Clearly M™ u is real symmetric by cyclicity of the trace 
and the fact that \JE m T v \/E m and t m are hermitian. (Ac- 
tually ip(A) is real for any complex d x d matrix A, and 
real symmetric for any A hermitian). Let v = (u ) £ M. d . 
Using ijTUjl we get 

V T ^{^fEm)v = V^M™V V = ^V^iV^TuVE^T^ 
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This follows from Lemma 0] since v^t^ £ Hermd(C). 
Hence M™ is a positive real (symmetric) matrix. 
The properties on purity simply follow from general facts 
on quantum operations on density matrices which remain 
true for generalized density matrices: 
For \u >< u\ and \v >< v\ generalized pure states, for 
any A complex d x d matrix and any generalized density 
matrix p, 

A\u >< u\A^ = \Au >< Au\ and 
\v >< v\p\v X v\ =< v\p\v > \v >< v\ (11) 

are generalized pure states. Relation follows from 
© and Tr(VEwV^) = Tr(E m p). □ 

The following properties will help to give a geometrical 
intuition of the action of the ip(s/Em) , s. For \/E m = 
\v >< v\ pure, < v\p\v >= Tr(V^Tp) = (l/d) ^/E^ .p. 
Thus using ifTTJl: 

— d 

So ip(\/E m ) is as was expected a non-normalized projec- 
tion. For any \JE m e HermJ(C) , the d 2 eigen- vectors 
v" of ip(^/E m ) with eigen- values X a correspond to 
d 2 hermitian matrices M a = 0~ 1 (wl) which satisfy 
V 'E m M° 'V 'E m = A a M cr (no summation). As a conse- 
quence, if p S C is such an eigenvector, then the rescaled 
density matrix p is such that p = p' m , i.e. p is unchanged 
if outcome m occurs. 

We now give the general expressions for p m in terms of 
p = (l/d)p M T M and y/E m = (l/d)yfe v T v , where we drop 
the index m and do not underline the components of the 
vectors p and \/E m for convenience. By definition: 



Pm = -p^fe^Pv^flcT^TvTc, 

Expanding this using tq = I and grouping the products 
of the Ti's in hermitian terms, we easily derive: 

Pro = V^oPoVeo 1 + ( 2 V£oPoV£i + V^oPiV^o)^ 

+ hV^iPoVEj + VeoPiVe^iTiTj + TjTi) 
+ ^V^tPjV^kinTjTk + TfcTjT;)} (12) 

To push the general d-dimensional analysis further, we 
need a particular choice of r^'s whose anti-commutation 
relations are convenient. This is subject to current work. 
We now treat in full detail the d = 2 (one qubit) case and 
apply our geometric approach to a challenging quantum 
information theoretical problem. 

III. THE QUBIT CASE PUSHED FURTHER 

Applied to qubit states the representation yields two 
of the most familiar objects in fundamental physics: the 



2x2 density matrices yield a Minkowskian future-light- 
cone in E 1,3 whose vertical sections are nothing but Bloch 
spheres. The correspondence between light-like vectors 
and fully determined spins is puzzling, but it requires 
a little more than one qubit to be investigated further. 
Meanwhile in this simple case we are able to give explicit 
coordinates for states posterior to non trace-preserving 
quantum operations. These formulae remain simple pro- 
vided Minkowskian products are introduced alongside 
the Euclideans. They constitute a sufficient armoury to 
deal, using only four-vectors, with the most general evo- 
lutions to happen on a qubit. 

A. The Cone and the Bloch Sphere 

A suitable Hilbert-Schmidt orthogonal basis for 2x2 
traceless hermitian matrices is given by the set of Pauli 
matrices: 

Tl = x= (? o) 

Together with the identity 

(J!) 

one may express any 2x2 hermitian matrix as a sum 
A = |j4 Tfj, with the ^4^'s real. The positivity conditions 
for those matrices turns out simple. 

Lemma 5 The cone of positive hermitian matrices 
Herrn^ (C) is isomorphic to the following cone of revo- 
lution in M 4 ; 

3 

r = {(A M )eM 4 /Ag-^A2>0,A >0} 

i=l 

Generalized pure states lie on the boundary ofT. 

Proof : The eigenvalues of A are given by A± = i (A ± 
y/AiAi)- A is positive if and only if A + A_ > and 
A+ + A_ > 0. This is equivalent to: 

> A A > (13) 

with ?7 M „ = Diag(l, — 1, —1, —1). The purity condition is 
an obvious consequence of Proposition ^ □ 

Thus the generalized (not necessarily normalized) den- 
sity matrices of a qubit cover the whole Minkowskian 
future-light-cone in E 1,3 . Taking a vertical cross-section 
of the cone is equivalent to fixing the trace Ao °f the den- 
sity matrix, which might be thought of physically as the 
overall probability of occurrence for the state. By doing 
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so we are left with only the spin degrees of freedom along 
X, Y, Z, therefore each vertical cross-section is a Bloch 
sphere with radius a — Aq. 

The ability to represent states with different traces 
is convenient when dealing with quantum ensembles 
{{px, Px)}x- When we seek to represent non trace- 
preserving quantum operations the feature becomes ab- 
solutely crucial. 

B. The Post-measurement State 

As we have seen in subsection III CI the most gen- 
eral quantum operation can be described as {M m } m = 
{U m y/E m } m with U m unitary and s/E m positive (the 
only exta feature Kraus operators allow is the possi- 
bility to ignore one's knwoledge of some measurement 
outcomes, but in our setting this is easily catered for 
by adding up the undistinguished non-normalized post- 
measurement states). While the action of U m is well 
understood in terms of four-vectors (as a mere rotation 
in the Bloch Sphere, see Lemma EJJ, the authors of this 
paper are not aware of a solid geometrical framework for 
representing the effects of y/E m - other than the one pre- 
sented here. In Lemma if A = \/E m while p is the 
initial state, then ApA stands for the (not renormalized) 
'post-measurement' state when outcome m has occurred 
(up to a unitary evolution U m ). 

Lemma 6 Let A and p be two matrices in Herm£(C). 
Then: 

ApA =i[-^(v44,,) + 2A {A.p)]T () 

+ -leaivw'AA*') + 2 A 3 (A-p)^ 

=l[r)»»>P v (r)w>A„A^) + 2A v ,{A.p)\r lJ , (14) 
Proof: Consider 

A= [a /3 7 6] 
p=[a x y z] 

We have: 

ApA =1 [a(a 2 + f3 2 + j 2 + 5 2 ) + 2a((3x + ■yy + 5z)]t 

o 

+ \[x{a 2 + [3 2 - 7 2 - 5 2 ) + 2(3{aa + 7 y + 5z)]n 

o 

+ \{y{a 2 - (5 2 + 7 2 - 5 2 ) + 2 1 {aa + f3x + 8z)]t 2 
8 

+ ^[z{a 2 -P 2 - 1 2 + S 2 ) + 2S(aa + f3x + ly )]r 3 

(15) 



This formula can be be obtained either by brute force 
calculation using the Pauli multiplication relations, or 
by exploiting the fact that Pauli matrices form a Clif- 
ford Algebra I.e. {rj, Tj} — 2SijTQ together with equation 
(|12[) . Regrouping the terms gives formula (|14fl . □ 

Corollary 1 Let A and p be two matrices in Herrn^ (C) . 
ApA can be expressed as a linear combination of p, A 
and the Identity: 

ApA = ^(A.p) A + ^{n^'A^A^ip - p T ) 

This last corollary provides much geometrical insight on 
non trace-preserving quantum operations. We find that 
the effect of \JE rn is not that difficult to visualize: the 
resulting state is a weighted sum of y/E m , the initial state 
and the identity, with real coefficients. 

It is a somewhat strange fact that the structure equa- 
tion i|15|) does not become apparent until one brings the 
Minkowskian product to the rescue. The spurious ap- 
pearance of special relativistic products in quantum me- 
chanics bears some explanation in this setting however, 
since the Minkowski metric is intrinsically related to the 
characteristic function of pointed cones of revolution. 

Finally it is important to notice that the results ex- 
pressed in these two last subsections are invariant under 
any orthogonal change of basis {r^ji. This is because 
rotations about the vertical axis leave the Minkowskian 
product invariant. The Pauli matrices have been help- 
ful in computing those results, but from now and in the 
rest of the paper we may consider ourselves in the more 
general setting of section HT1 

C. Square and Square Root 

In our quest towards representing non trace-preserving 
quantum operations in the cone we have managed to ob- 
tain the probability of occurrence p(m) in terms of E_ m 
(Proposition^. In the previous subsection we have also 
worked out the evolved state p m , but unfortunately this 
was done in terms of y/E m . In order to deal fully with 
these operations in the Cone formalism we need to un- 
derstand ways of switching back and forth from E_ m to 
y/E m . The next Lemma is a direct consequence of equa- 
tion (| 1 4|) when p = I. 

Lemma 7 The square of a matrix A in Herrn£(C) is 
given by: 

A 2 = AoA-^iv^A^To 
Inversely the square root operation obeys: 

VI = -(A+ -yti^ro) 



with: r = 
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Note that A is proportional to V A if and only if A is 
generalized pure or A oc I. 

But when we seek to express a function of E_ m in terms 
of \JE m (or the reverse) the next formulae become con- 
venient. 

Lemma 8 Let A and p be two matrices in Herm\ (C) . 
The following relations hold: 



Proof: By using @ we have: 



A 2 -p = MA-p) - 2fi)(v4(A) 



with: r = ^Aq + ^n^vA^Av 



On the whole taking the square root of E_ m is not so easy. 
It would be much more convenient if we could make all 
calculations in terms of E_ m , with the added advantage 
condition: 



(16) 



is easily visualized. Results in the following subsection 
are most useful for this purpose. 



D. Inner Products Through Quantum Operations 

Consider two states p°, p 1 . Suppose they undergo a 
quantum operation {M m } m — {U m \/E m } m and outcome 
m occurs. Rather than seeking the coordinates of the 
rescaled post-measurement states p® m and p Y m , we may 
be interested in their positions relative to one another. 
Note this subsection reuses a number of notational con- 
veniences introduced in section ITT1 

Lemma 9 Let p° ,p x be two initial states in Herrn^ (C) 
I Em a measurement element in Herrri£(C). The 



inner products of the post-measurement states satisfy: 



P° P 1 

rm rm 



(17) 



' 1 ' 
Pm 'Pm 



P° P 1 
rm rm 



vv' P V P v , ) 
(Ern.pP)(Ern.pl) 



-[(En.ftiEn.f) - (r^^^H^p^,)] 



' 1 
Pm 'Pm 



(E^.p^iE^.pi) [ 1 



Pm ' Pm 



= 2Tr(^E~p°^E~^EZp 1 VE^) 
= 2Tr(E m p a E m p 1 ) 
= E m p°E m . p 1 

From there we readily obtain equation (|17l) by applying 
equation l|14|) once. □ 

By letting p° = p 1 = p in the above Lemma we get: 

llPmll 2 = \\2{E m .p) 2 - (r)^E m ^E mfi , ) (j] vv >p v p v , )] 

(l W .'£m,£m„,)(V/)iJ 



IKjf 

ll/^ll 2 



{Ern.p) 2 

^[(Ern.p) 2 - (r) uu , E m ^ E m/j , ) (r\ vvl p r p r , ) 

{Era.pf 



(19) 



Equation (|19|l clearly exhibits the general property we 
stated in Proposition [21 that is if the initial state is 
generalized pure (^^'P p , = 0) or the measurement 
is generalized pure ( r) uu ' E m E m , = 0) then we have 

IIPm'H = 1 (pure), which implies that p m is generalized 
pure. 

The above lemma enables us to determine all the rela- 
tive positions (angles and norms) of quantum states using 
relatively compact formulae which do not involve \J E rn . 
It is only when the coordinates of each post-measurement 
state are required that one needs to take the impracti- 
cal square root of E m . But remember we are allowed 
an arbitrary rotation U m in order to complete the quan- 
tum operation. This means we have full freedom to fix 
the absolute coordinates at will (so long as the relative 
positions are respected). 

Most Quantum Information Theoretical problems seek 
to evaluate the limits of quantum operations, e.g. quan- 
tum cloning distinguishability |10| . Infomation Gain 
versus Disturbance tradeoff In these situations the 
precise individual coordinates of the states after Ad^^ 
tend not to matter; usually they will need to be rotated 
anyhow into a position which optimizes the fidelity mea- 
sure in question. What counts is the relative position of 
the post-measurement states. Therefore these problems 
can be treated comfortably in our framework. Section HVI 
provides a good example of such an application. 

There are, however, some rare situations where we 
would like to see quantum operations act step by step, 
yielding precise coordinates - instead of just fixing the 
coordinates of the final state as we would do in order to 
avoid taking the square root of E m . This is the case for 
instance in quantum complexity, where one needs an ap- 
preciation of how many basic computational operations 
it takes to accomplish some calculation. Yet in this type 
of problems it turns out that the basic operations can be 
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taken to be unitary operators, with measurements only 
performed at the end (principle of delayed measurement 
Therefore these scenarios may still be analyzed com- 
fortably within our conal representation: the basic uni- 
tary operators will just be a set of chosen real orthogonal 
rotations, and the final measurement statistics will be 
evaluated straight from E m . 



A. Information Contribution, Disturbance 
Contribution 

Suppose the {\ipx)} x =o,i states Alice prepares verify 
the following basic relations: 

Vl= <t>(\*p x )& x \) 
l£_.vi = d= \J\ - c 2 



By choosing a suitable basis in the Bloch Sphere and 
since the are pure we may fix: 



IV. APPLICATION: INFORMATION GAIN 
VERSUS DISTURBANCE TRADEOFF 



The following idealized scenario captures a key situa- 
tion for any quantum cryptographic protocol: 

Alice owns a random variable X = {(|,0),(i,l)}. Ac- 
cording to the outcome x she prepares either l^o) or 

i.e. she runs |ir)|0) \x)\ip x ). Eve knows U and the dis- 
tribution X, but not the particular outcome Alice has 
drawn. Later Eve gains access to \tp x ) and may use of 
this opportunity to try and learn about x. How much 
she learns is quantified using Information theoretical no- 
tions. Even though Alice has had to expose \ip x ), still 
she really wanted to keep x secret. But now she gets a 
chance of checking upon Eve's honesty - by asking her 
to return \ip x )- Suppose Eve's measurement and further 
manipulations have modified \ip x ) into p x . Alice then 
measures {\ip x )(ip x \^ Id — |'0a;)('0x|} and has a probabil- 
ity 1 — (ijj x \p x \ip x } of detecting the felony. 

The point is that most quantum cryptographic pro- 
tocols rely upon the fact that Eve cannot eavesdrop a 
state without causing it an irreversible, detectable dam- 
age. In spite of their central role, Information Gain ver- 
sus Disturbance tradeoffs upon discrete ensemble states 
remain largely unknown, due to the mathematical diffi- 
culties they raise. In 1995 Fuchs and Peres accomplished 
the mathematical feat of obtaining an analytic formula 
for the above case of two non-orthogonal states. But the 
method they used relies upon a number of "plausible" 
assumptions - and does not provide a geometrical intu- 
ition of what the family of optimal measurements looks 
like. 

The Cone, by enabling a per outcome geometrical rep- 
resentation of generalized measurements, permits us to 
overcome some of these shortcomings and greatly facili- 
tate the derivation of Fuchs and Peres' formula. We hope 
this illustrates the power of the geometrical framework 
developed in this paper. 



rP = [ 1 c dO 
yi=[l -c d 



(20) 
(21) 



The most general thing Eve can ever do is to attack 
the states with a measurement {M m } m . This procedure 
is equivalent to first measuring {\/E m } m , and then, con- 
ditional to to, applying the unitary transformation U m , 
with E m and U m defined as in subsection III CI 
It is rather interesting to observe that the second step 
has no other use but to "repair" the post-measurement 
states as much as is possible. The first step on the other 
hand may partially destroy the initial states so as to col- 
lect the Information Eve seeks. This is the step we now 
study in order to quantify her Information Gain. 

Let Y be the random variable arising from the mea- 
surement outcomes, i.e. Y = {(p(m),m)} m . Eve's Infor- 
mation Gain is given by: 

/ = H{X : Y) = H{Y) - H{Y\X) 
— ^]p( TO ) l°g(f>(m)) — y^p(z, to) \og{p{m\x)) 

m X ,m 

= I m with 

m 

I m = p(m) \og(p(m)) - ^p(x, to) log(p(TO|x)) 

X 

I m must be understood as the Information Contribution 
brought by the measurement element: 

= [ a (3 7 5 ] = (f>(E, n ) 

By making use of the relations H2()|l , H21JI and one can 
express I m geometrically in terms of scalar products in 
the cone: 

I m =- (Pm + Qm) log(p m + q m ) 
log(2p 

in ) + g m log(2g m ) 

a + PC + id Ern-V^ _ 

with p m = = = A — —p(v,m) (22) 



4 4 
a - f3c + ~fd e m -v^_ 



= p{l,m) (23) 



Notice that if Era is orthogonal to v}_ (resp. v^_) then 
Im = Pm (resp. q m ). Such a measurement element may 
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be said to be "all or nothing" : it brings a whole bit 
of information when it occurs, but does so only with 
probability p m (resp. q m )- Taken individually these 
measurement elements seem ideal: they fully identify 
\ip x ) and thus they let you reconstruct the initial state 
perfectly, with no disturbance at all. The downside is 
that failure to occur comes at a high price. In order 
to verify the condition the other measurement 

elements generally become rather inefficient with respect 
to the tradeoff. The family of the optimizing {M m } m is 
not constructed in such simple ways. 

Next we seek an expression of the Disturbance Con- 
tribution brought by each measurement element. For 
this purpose we must first assume outcome m has oc- 
curred. Eve knows it, and now she will try to maximize 
her chances of fooling Alice by applying a carefully tai- 
lored unitary evolution U m . First we will give D m as a 
function of U m , and next proceed to the maximization 
which determines U m . Remember that upper indices x 
distinguish initial states, while lower indices to specify 
the measurement outcome. 

p(fool|m) = Y^P{x\m)Tz{\^ x ){il) x \U m p x JU m ) 

X 

T,xP( x \ m )^- r m 



of U m we can safely assume rj^ , r m , v° , v 1 to be coplanar. 
Thus D m may now be rewritten in terms of those angles 
as well as p m and q m : 



Pn 



p m COs(A m - LU rn ) - q m COs(A m + LU m ) 



D, 



In this equation the values of p m , qm and A m are fully 
determined by e m , as described in (|22[) . (|23|) . H17[) . cu m on 
the other hand solely depends on U m : it can be chosen 
at will by rotation in the Bloch Sphere. We now show 
how Eve must tune iu m so as to minimize D„, . 



dDn 

duj„ 



=> p m sin(A m - Lu m ) - q m sin(A m + to m ) = 



The minimum occurs at: 



arcsm 



' \Jv m + <iL + 2Pmq m cos(2A m ) 
which yields, after simplification: 



Pm T" Qm 

Vp m + i m + 2 Pmq m cos(2A m ) 



B. The Tradeoff 



1 + T lx p( x \ m ) vX - r m 



where 



L m = [l rt ] = m m p X m'Um) 

_ (j>{U m y/^\^x){^x\y/EmU m ) 
p(m\x) 

Negating back to the Disturbance we obtain: 
D = ^2 D m with 

m 

D m = p( _, fool, m) 

_ P(m) - Y, x P( x ' m^.r^ 



= P( m ) - J2 x P( x ' m )\\ vX \\ \\ r rn\\ C0S (^:O 
2 

In our scenario the {|"03;)}a: are pure. Thus by Lemma|21 
or equation (|19fl we have ||r^J| = 1. Now let us deal 

with cos (v x , r m ) by making the following definitions: 

6 = $y) 



How many elements should Eve's measurement con- 
tain? Lcvitin has proved that there exists a two-element 
measurement jM m }m=o,i which maximizes Eve's Infor- 
mation Gain [lj. While this was never formally shown 
to be the case for the measurements which optimize the 
Information Gain versus Disturbance Tradeoff, there is 
strong numerical evidence in support of this assumption 
@. Suppose this is the case and let e m denote the jjL th 



coordinate of 
have: 



V/ 



Using the constraint equation (|16fl we 



Se . 



-5s i and thus: 



df df 



de , 



— m 



(24) 



Optimizing the Tradeoff implies finding a stationary 
point for the Disturbance while keeping the Information 
Gain fixed. We need to find £o such that 



^^ £0 =0 



^ de —>• 
n — n 



@m — { r rrn r m) 
A m = 8 — 9 rn 

Lu, n is the angle between the bisector of (r m , rj n ) and that 
of (y , v 1 ). Given that we want to minimize D m in terms 



where the variations feo M are subject to the additional 
constraint: 

|U — M 
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Using equation (|24() and D = Dq + D\ and I = Iq + Ii 
this gives: 



subject to 



/J — M /J — M 

\ - dl0 r 

|U — — M ju —A* 



(25) 



(26) 



Guided by the geometrical picture of the scenario one 
may consider the following attack: 

£0 = [1 P 0] 
ei = [l -/3 0] 

The fact that this is indeed a solution follows from its 
obvious symmetries: 



For n ^ 2 



For /i ^ 2 



<9eo der 



— and tor p = 2 



9eo 9ei 



and for p = 2 



9eo, 



(27) 
(28) 



may undergo is a generalised measurement {M rn } rn = 
{U m -\/E m } m , where the polar decomposition was ap- 
plied. In order to represent M m 's per-outcome effect 
upon the real vectors we defined ip : A <fr o Ad a o^ 1 
and showed that tp(U m ) is a real orthogonal transform 
while ipiVEm) turns out to be a real positive matrix. 
Thus the geometrical effect of a generalized measurement 
can be viewed in terms of real transformations only. 
Such a nice correspondence suggests quantum mechan- 
ics could be expressed elegantly over the real numbers 
in this manner, quite differently from its formulation in 
terms of real Jordan algebras 12|. However we first 
need to gain more geometrical intuition about the set 
of real vectors <^(HermJ(C)), and the sets of allowed 
orthogonal and positive transforms. For now we know 
that </>(HermJ(C)) is a subcone of the future-light-cone 

r = {(a m ) e R^/zSi 1 a? < (d- i)Ag, a > o}. 

One of the advantages of defining <\> upon HermJ(C) 
instead of the restricted set of density matices is that 
E m = MmMh can be visualized. In order to character- 
ize its effects we derived rather compact and powerful 
formulae for the qubit case, such as the one giving the 
scalar product of the post-measurement states: 



Substituting 128|) in the constant Information constraint 
we get feo 2 = 0. Using this fact together with equa- 
tion H27fl it becomes clear that condition l|25[) is fulfilled. 
Thus £0 is a stationary point. We may now proceed to 
compute the values of the Disturbance and the Informa- 
tion Gain under this family of optimal attacks. First by 
making a few additional observations: 

Po = Qi = P 

Pi = qo = q 

D = D 1= D/2 
I =h = 1/2 
1 



D = |- y/p 2 + q 2 + 2pq + cos(2A m ) 
J=l + 2plog(2p) + 2glGg(2g) 

and second by plugging in the relations l]18p. (|20|) - (|23|l . 

we reproduce the exact content of Fuchs and Peres' for- 
mulae: 



D = ~ - ly/l + (c 2 - c 4 )(/3 2 - 2 + 2^1^) 
I = + /3c) log(l + /3c) + (1 - fie) log(l - /3c)) 



V. CONCLUSION 

In this paper we considered a linear embedding tak- 
ing d x d positive hermitian matrices into vectors of d 2 
real entries, <j> : p 1— > p — (Tr^r^))^. It is a well-known 
fact that the most general evolution a density matrix p 



By looking at such expressions it becomes apparent that 
Minkowskian products have a crucial role to play in our 
framework, and even more so as we showed that pure 
quantum states correspond to light-like vectors (i.e. they 
sit on the boundary of T), even in dimensions greater 
than 2. It seems interesting to notice that BPS states 
of supersymmetric theories can also be thought of as ly- 
ing on the boundary of a cone of positive operators |?|], 
and that their stability is related to that fact. Somehow 
our setting seems to single out generalized pure states in 
a more natural way than merely characterizing them as 
unit rank elements of the boundary of Herm J (C) . The 
reminiscence of special relativity must be investigated 
further; this will be a subject for future work. 
Pauli matrices together with special relativistic consider- 
ations have already brought some fruitful results to quan- 
tum information theory. This is the case for instance in 
where some limits of quantum cloning are derived 
by using the no-signalling condition. Armed with the 
present representation one should be able tackle more of 
these difficult quantum information theoretical problems. 
Already in this paper we recovered Fuchs and Peres' in- 
formation gain versus disturbance formula simply and 
geometrically. In the future we should be able to extend 
our analysis to the case of two non-equiprobable states. 
Some highly symmetric n-states scenarios may well cease 
to be out of reach. 
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